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ABSTRACT 

Using  electromagnetic  waves  to  detect  and  identify 
complex  structures  near  the  earth’s  surface  remains 
challenging  because  of  ground  conductivity  and  ground 
clutter.  Powerful  computer-based  numerical  techniques 
are  required  to  solve  canonical  propagation  problems 
involving  electric  dipoles,  magnetic  current  loops, 
perfectly  conducting  and  dielectric  spheres,  cylinders,  etc. 
A  critical  need  remains  for  real-time  realistic  and  robust 
analytical  techniques  that  extend  detection  and 
identification  domains  to  complex  structures.  Our  study 
addresses  this  need  using  electromagnetic  reciprocity  in 
combination  with  dyadic  Green’s  functions  and 
asymptotic  theory  of  integration. 

1.  INTRODUCTION 

In  the  past  decade  theoretical  interest  for  detecting 
and  identifying  surface  and  shallow  buried  complex 
targets  has  been  minimal  to  modest  in  favor  of  more 
expedient  experimental  studies.  As  recently  reported  at 
the  international  EUROEM  conference  (Lausanne, 
Switzerland  July  2008)  success  in  meeting  desired  useful 
operating  curves  for  mine  detection  has  not  been  entirely 
successful  using  various  combined  theoretical  and 
experimental  schemes  for  determining  the  Probability-of- 
Detection,  Pd  ,  and  Probability-of-False  Alarm,  Pfa  . 

Two  reasons  given  for  the  slow  lack  of  progress  in 
determining  Pd  and  Pfa  are:  (1)  this  is  a  fundamentally 
difficult  detection  problem  with  true  "natures  limits”  on 
Pd  and  Pfa ,  and  (2)  we  have  pushed  our  theoretical 
capability  to  its  limits.  It  is  our  perception  rather  that: 
(1)  we  haven’t  even  come  close  to  exhausting  detection 


algorithms,  and  (2)  the  powerful  technique  of  dyadic 
Green’s  functions  (Tai,  1994)  has  not  yet  been  applied  to 
detection  near  the  ground-air  interface.  In  this  paper  we 
address  the  latter  idea.  In  Section  2  we  define  the  problem 
to  be  solved.  The  theoretical  model  based  on  the  dyadic 
Green’s  function  approach  is  described  in  Section  3. 
Detection  theory  using  the  integrated  concepts  of 
reciprocity  (Monteath,  1973)  and  the  asymptotic  method 
of  integration  (Chew,  1990)  are  rendered  in  Section  4.  A 
calculational  approach  and  sample  problem  are  discussed 
in  Section  5. 


2.  STATEMENT  OF  PROBLEM  AND  OVERVIEW 

Figure  1  depicts  the  detection  and  identification  of  a 
buried  object.  In  a  related  paper  von  Laven  et  al.  (2008) 
have  presented  experimental  results  dealing  with  this 
subject.  The  physical  basis  for  doing  this  began  nearly 
one  hundred  years  ago  with  the  observation  that  the 
ground  played  a  major  role  in  the  propagation  of  radio 
waves.  Arnold  Sommerfeld  (1909)  provided  an  elegant 
theoretical  model  of  how  ground  conductivity  affects  the 
propagation  equations.  These  equations  turned  out  to  be 
so  challenging  that  special  mathematical  approximations 
were  developed  for  the  case  when  antennas  are  located  on 
the  earth’s  surface.  In  the  mid-1930s  Norton  (1937) 
developed  robust  and  yet  practical  models  of 
electromagnetic  radiation  from  horizontal  and  vertical 
electric  and  magnetic  dipoles  near  the  earth’s  surface. 

During  WWII  and  moving  into  the  early  1970s 
research  progressed  on  buried  antennas  of  various  sorts 
and  on  underwater  communications.  Ground  penetrating 
radar  (GPR)  surfaced  in  the  1970s  when  utilities  tired  of 
continuously  repairing  underground  water  lines,  gas  lines, 
and  electric  power  cables.  GPR  is  generally  successful  in 
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finding  these  relatively  large  targets.  In  the  early  1990s 
and  continuing  today  there  is  great  concern  for  detecting 
large  buried  objects  such  as  bunkers  deep  into  the  ground, 
but  this  is  different  in  quantitative  ways  from  detecting 
the  surface  and  shallow  buried  objects  of  interest  in  this 
study. 


Figure  1.  Detecting  and  identifying  buried  object. 


Our  targets  differ  significantly  from  canonical 
electric  and  magnetic  dipoles;  detection  and  identification 
at  suitable  range  is  also  critical.  Kraichman  (1970)  has 
presented  a  comprehensive  treatment  of  buried  dipoles, 
and  Iordan  &  Balmain  (1968)  have  presented  a  simplified 
version  of  propagation  models  that  include  practical 
curves  of  attenuation  over  ground. 


ground.  Using  Tai’s  dyadic  Green’s  function  notation,  the 
relationships  between  the  electric  fields  and  currents  are 

E2(f)  =  ioj^  j  Gm>  (r2  7  >  7,  (f '  )dVt  1  (2a) 

V/ 

Eff^^ioj^  J  G(12)(FJF2').72(F2W2'  (2b) 

V2‘  1 

The  dyadic  Green’s  function  Gf2i)(r2\r\)  connects  the 
electric  field  below  the  earth’s  surface  to  the  source  above 
the  earth’s  surface  while  Ga2l(F,|F2')  connects  the  field 

above  the  earth’s  surface  to  the  ground  current.  The  “  •  ” 
symbol  is  the  usual  dot  product.  Depending  on  the 
properties  of  72(F2 ')  it  may  be  more  efficient  to  transform 
the  right-hand-side  of  Equation  (2b)  into  the  following 

form  derived  from  the  reciprocity  principle  (the  sign 

means  transpose). 

G(12,(F1|F2,)=G<21,(F2’|F1)  (3) 

g,12,(f1|f2,).72(?2')=g(21,(H)'|f1).72(h)')  = 

1  L  (4) 

72(F2').G,21,(F2'|F1) 

The  foregoing  reciprocity  results  are  easily  derived  for 
our  case  where  //,  =  /J2  =  jU0 . 


The  detection  process  comprises  the  following  steps: 
(1)  generating  the  below-ground  electromagnetic  field 
from  the  transmitter,  (2)  generating  a  below-ground 
current  via  target  interaction,  and  (3)  determining  the 
aboveground  electric  field  from  the  below-ground  current. 
In  most  cases  of  practical  interest,  estimating  the  below¬ 
ground  field  from  an  aboveground  transmitter  operating 
in  the  100  MHz  and  above  frequency  range  is  relatively 
easy,  requiring  mainly  knowledge  of  the  ground  reflection 
coefficients.  Predicting  the  generation  of  intrinsic  target 
current  and  related  target  induced  ground  current  is 
difficult  because  detailed  knowledge  of  targets  is  not 
always  available. 

We  select  a  generic  target  class  characterized  by  the 
constraint:  jjjJ2(r')dV'  =  0 .  This  equation  means  that  the 
current  is  generated  wholly  within  the  target  with  each 
component  satisfying  the  equation 

tflJ2«(f')dV'=0  a  =  x,y,z  (1) 


A  major  building  block  for  solving  the  two-region 
system  of  Fig.  1  is  to  use  the  free  space  dyadic  Green’s 
function  given  by 

G0(F  I F')  =  fj  +  (l/k2)w)G0(F  I F’)  (5) 

In  Equation  (5)  the  dyadic,  7 ,  is  the  idem  factor 
( ax , ay , dz  are  unit  vectors),  and  G0(FIF')is  the  free 
space  Green’s  function. 

I  =  axax  +  ciyciy  +  a.az  (6) 

G0(F  I F')  =  exp/k  I F  —  F'l / 4zr I  F-F'l  (7) 

In  summary,  this  paper  shows  how  to  derive  practical 
expressions  for  the  aboveground  field  from  target 
currents.  Equation  (2b),  from  the  mathematical  properties 
of  the  dyadic  Green’s  function  and  the  asymptotic 
approximations  to  complex  integrals. 


We  use  the  following  notation:  a  prime  symbol,  such 
as  r\  refers  to  a  source  coordinate  while  the  unprimed 
coordinate,  such  as  r,  refers  to  a  location  at  which  the 
electromagnetic  field  is  computed.  The  subscripts  1  and  2 
refer,  respectively,  to  sources  located  in  the  air  and  in  the 


3.  SUMMARY  AND  REVIEW  OF  DYADIC 
GREEN’S  FUNCTION  APPROACH 
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In  this  section  we  summarize  the  dyadic  Green’s 
approach  for  determining  the  electromagnetic  field  from 


sources  located  above  and  below  the  ground,  respectively. 
This  material  is  extracted  from  the  excellent  text  by  Tai 
and  is  provided  to  allow  the  reader  to  apply  the  method  to 
related  scenarios.  Determining  the  electromagnetic  fields 
for  Fig.  1  is  best  achieved  in  cylindrical  coordinates. 

The  starting  point  is  to  determine  the  free-space 
dyadic  in  cylindrical  geometry  from  the  equation 

VxVxG0(FI  r')-k2G0(r  \r')  =  IS(r  -r')  (8) 


where  IS(r  -  r' )  is  expressed  in  the  following  orthogonal 
set  of  basis  cylindrical  wave  vector  functions: 


MenA  = 


MonZ  = 


nJJAr)  .  dJ(Ar) 

- - sin  n(par - - - cos  ncpa 

r  dr  ' 


,  nJ„(Ar )  dJ  (A  r)  . 

H - - - cos  n<par - — - sin  ntpa. 


Ne,a  k 


N o,a  ~  , 

*■1 


.,  8J n  (A  r)  ihn  ,  . 

ill — - cosnipar - J  n(Ar)smn<pa, 

dr  r 

+  A2Jll(Ar)cosn0az 


8J „(Ar)  .  ihn 

ill - - - sin ntpar  H - J  n  (Ar)cosntpa , 

dr  r 


+  A2  Jn(A  r)sinn0az 


kj=A2+h2 


exp  (ihz) 
(9a) 
exp  (ihz) 
(9b) 

exp  (ihz) 

(9c) 

exph'fe) 

(9d) 

(9e) 


In  the  foregoing  equations  J  n  (At')  is  an  integer  Bessel 
function  of  order  n. 

Equations  for  the  two-region  problem  of  Fig.  1  are 


VxVxEl-k2El  =i<DJu0Jl  z>  0  (10) 

Vx  VxE2 -k2E2  =0  z<0  (11) 

At  the  ground-air  interface  we  have 
HxE1=nxE2  (12) 

nxHl=  n  xH,  (13) 

VxE=ico/j0H  (14) 

V  xH  =-ia>sE  (15) 


nx'VxH,  nx'VxHn 

- : - = - : -  (16) 

-icoel  -i(oe2 

nxVx£j=nxV  xE2  (17) 

When  the  source  lies  above  the  ground,  we  use  the 


following  Green’s  functions: 

VxVxG*111  -k2G2ll)  =IS(r-r')  z>  0  (18) 

VxVxGf"  ~k2G(21)  =0  z  >0  (19) 

If  the  source  is  in  the  lower  region,  we  have 

VxVxGf-  k2G212)  =0  z  <  0  (20) 

V x Vx Gj22)  -k\G(2T)  =  I5(r  —  r')  z< 0  (21) 

For  the  z  >  0  case  we  have 

Efi)  =  io^\Gfl\7  \?')»Jl(r')dVl'  (22) 

v,' 

£,(?)=  iwjio  JG<21)  (r  I r')»Jl(r')dVl'  (23) 

V,' 

For  the  z  <  0  case  we  have 

£,(r  )  =  icon o  J  G\l2)(r  I r  ')  •  J2(r')dV2'  (24) 

V2’ 

E2(r  )  =  iwjUo  |  Gf2)(r  |r  ')•  J2(r  ")dV2'  (25) 

VA' 


From  the  relationship  D  •  C  =  C  •  D  we  get  the 
equivalent  ways  of  computing  E{(r  )  and  E2(r  ): 


E,(F)  =  i&tt,  J  ')  •G^11)(?  I r  ')dV  '  (26) 

Efi  )  =  toft  J  7,(7  ')  •G<1I)(r  '\r  )dV  '  (27) 

vz 

E2(r  )  =  ico/u0  J  7, (7  ')  •  Gj21)(7  |r  ')dV  '  (28) 

v2' 

E2(r)=  ia>, Ho  J  7,(7 ')  •  G? '>(7  '\r  )dV  '  (29) 
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The  building  blocks  for  the  solution  of  the  different 
Green’s  functions  are  the  wave  vector  functions  used  in 

conjunction  with  G0(r  IF'). 

In  region  1  with  z,z'>  0  and  z>z\  we  have: 


Z-  /  oo  oo  9  _  A 

G0Cr\r')  =  —  ldAZ^jTJL 

47T  o  n= o  An, 


MenAh1)M'e„A~hl)  + 

MolU(hl)M'oia(-h1)  + 
NenAhlWmA~h1)  + 
N0MN'onA-h ,) 


(30) 


In  region  1  with  z,z'>  0  and  z,’>  z  ,  we  have 


-  /  oo  oo 

G0(r\n  =  ^-!dAX—JL 

47T  o  n= o  An, 


MenA-K)M'mAh  i)  + 
MnnA(-hl)M'onA(hl)  + 
NenA-KWenAK)+ 
NonA~KWonAh ,) 


(31) 


aMmA\ W'mA  0h)  + 
aMonA(hl)M'onA(hl)  + 
bNmAK )N'enX  (/;,)  + 
bN  onAKW  onAK) 

(39) 

C^erA~h 2W'mX  (/!,)  + 

cMonA~h  2)M'onA(hl)  + 
dN  enA~h2W  CTl  (/»,)  + 
du  onA-hiW  onAK) 

(40) 

The  unknown  constants:  a,  /?,  c,  and  c/  are  found  by 
matching  the  boundary  conditions  at  the  interface,  z  =  0 . 
We  have 

nxGj111  =FxGf'  (41) 

iixVx  Gj1 11  =  //  x  V  x  Gf  ’  (42) 


-  /  CO  OD  7  —  A 

GL21)(?IF')  =  -^j^S-A 
4^  0  «=0  /i/lj 


g£1>(FIF')  =  —  JrfAI 


2-<5n 


4^r  o  n=o  T/zr 


The  prime  notation  used  in  Equations  (30)  and  (31), 
such  as  M'enl ,  means  that  it  is  a  function  of  the  primed 
coordinates,  r'  ,</>'  ,z'  ■  The  following  definitions  apply: 


c>0  =  1  ;  for  m  or  n  =0 

(32a) 

S0  =  0  ;  for  m  or  n  ^  0 

(32b) 

II 

S 

ii 

(33) 

k2  =  a>yjju0s(l  +  ia!  cos) 

(34) 

a  = 


lil  —hn 
h{  +lu 


A  hi  —  k2h2 
k2hi  +  k2h2 


2  hx 

h{  +  h 2 


2kjk2/ij 

k2hx  +  kfh2 


(43) 


(44) 

(45) 

(46) 


hi  =yjkf  -A2 

(35) 

J5" 

II 

tO  (o 

1 

(36) 

To  solve  Equations  (18)  and  (19)  we  assume 


4.  ELECTROMAGNETIC  FIELD  DETECTION 

The  model  described  by  Equation  (1)  assumes  no  net 
current  flowing  into  or  out  of  the  source  volume,  V.  All 
current  circulates  within  the  volume.  For  brevity  we 
assume  no  current  flow  in  the  z  direction.  We  then  have 


G‘11)(FIF')=G0(FIF')  +  G”1,(FIF')  z>0  (37) 

Gf1,(FlF')  =  G‘;1,(FIF’)  z<0  (38) 

In  order  to  satisfy  the  boundary  conditions  at  the 
interface  between  the  two  regions  and  the  radiation 
conditions  at  z  =  ±°°,  we  must  have 


j\J2x(x',y')dx'dy'=0  (47) 

\\J2y(x'  ,y')dx'dy'  =  0  (48) 

Even  though  the  foregoing  equations  are  satisfied,  we  still 
have  Jfx',y')  ^0,  J  v(x',  y')  ^0  for  f./eV  .  Using  the 
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relationships  x'=r'cosi/>'  and  y'=r'smif)'  we  write  the 
radial  and  angular  components  of  current  density  as 

J2r  O' ,  <jl )  =  cos  </>'  J2x  (r' ,  <// )  +  sin  <//  J  2y  (r' ,  </>’ )  (49) 

J  2^(r',  </>')=  COS  <fl'J  2.M'  +  J2y(r'  ’</>')  (5°) 

Using  the  unit  vectors  the  vector  current  density  is 

J2AJ)  =J2M\<tnar  +  /2p(r',^)a„  (51) 

The  aboveground  field  is  computed  from  the  source 
term  below  ground. 

£,(?,)  =  icm,  J  G<12)  (r, I r2 ')  •  72 (r2 ')rfV2 '  (52) 

V2' 

It  is  not  necessary  to  determine  GjI2)  (ij  IF,')  by 
repeating  the  calculations  of  the  previous  section.  All  that 
is  necessary  is  to  interchange  the  roles  of  ft,  and  ft, ,  and 

the  vertical  positions  of  the  transmitter  and  target,  Z1,Z2, 
respectively.  Thus,  all  the  previous  formulas  are 
applicable  with  the  following  substitutions  (new  values 
are  indicated  with  the  “hat”  symbol  such  as  a  ): 

ft,  =  co^/u0en  =  ft,  (53) 


With  this  upside-down  frame  of  reference,  we  replace 

G312,(rlF')  by  G^'VF  I  F') ,  and  we  let  r,(/>  denote  the 
field  point  of  interest  and  r\tji  the  source  coordinates.  To 
minimize  and  simplify  the  notation,  it  is  now  stated  that 
the  field  point  is  above  the  ground  (we  drop  the  subscript 
“1”)  and  the  current  is  the  target  current  (we  drop  the 
subscript  “2”).  We  still  retain  the  subscripts  for  the 
parameters  defined  by  Equations  (53)  to  (62).  Since  the 
height  of  the  true  field  point  is  the  positive  z, ,  it  is  now 
the  below-surface  fictitious  depth,  -z, ,  and  positive  z2  is 
now  the  height  of  the  source.  With  this  new  notation 
Equation  (52)  now  reads 

E(r  )  =  icoii0\Gfl\?  I?  '  W(F  ')dV'  (63) 

v  ' 

We  have 


cMmA(-h2  )M'enX(hl)  + 

ico^  =-](alC0^{2~5°) 

0  -0  A 7i Ah l  dN mk (~h2  )N'enA  (/z, )  + 

dN  ora(-h2)N' 


ft,  =  a>J  /u0s(l  +  ia  /  coe)  =ft. 


h2  =  yft,2  -  A2  =  /z, 


ft,  =  yft,2  -  T2  =  /z. 


ico^G?"  =-jdAZ 


cM„n2(-hl)M'enZ(h2)  + 
“  co^0(2-S0)  cM onX(-h{)M' onX  (/z,)  + 


o  n=o  4nAh2  dNenA{.-hl)N'era{h2)  + 
dN  0A~hiW  olU(h2) 


.  h,  -  h? 
a  =  -7T — 7T-  =  -a 


io)jU(iG<32^ (r  I F ')  =  io)^  £  G„ 


:  ft2/z,  -k\h2 
b  =  — — ^ .  =  —b 
k2h{  +  k[h2 


E(r)=  X  j Ni)jU0Gn  •  J (r')dV' 


2/z,  2/z, 

/z,  +  /z,  ft,  "t"  ft  i 

2ft,ft,/z,  2ft,ft,/z, 

ft,2  /z,  +  kfh2  kfh2  +  k2h{ 


E(r  )  =  S  En(r) 


imS.=4 

0  A- 71  An  ^ 


yVAn+%in+'VCn+^r>n 


An  =  cM enX(-h{)M' enX  (h2) 
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(75c) 


=cMonZ(-hl)M'on;i(h2) 


(69b) 


r^=-ja/,,ar,)cos  nfJ,dV' 
v  or 


'Vcn=dMm,{.-hl)M'enX{h2)  (69c) 


KVr>n=dMon/(-lh)M'aia(h2)  (69d) 


En  (?  )  =  Em  (? )  +  EBn  (r)+ECn(r)  +  EDn  (?)  (70) 


EAn  (?)  =  -f^(2„  ^o)  J**,  (?  I  ?'  W (?  W  (71a) 

0  4^>i/22  V’ 


ErJr)  =  -r^(2/o)^jfa„(r  I  r').J(r')dV'  (71b) 
o  47r/t/29  v 


4,(?)  =  -J^f„/^?iTc,.(F  I  r')*J(r')dV'  (71c) 
o  4;ta/29  v 


(?)  =  - J  {\  ,S° )  dl  jyD„  (?  I  ?' )  . ./(? '  W  (7 Id) 

o  4;r/t/i2  v 

By  example,  we  carry  out  the  calculation  for 

EAn(r). 


Since  there  is  no  vertical  component  of  volume  to 
consider  (in  this  case  J and  J  r.  have  the  dimensions  of 
amps  per  unit  meter),  we  have 


dV'=  r'  dr'  dip' 


U)  Lr2?nJMr')~. 


C=-JJ 

oo  r 


sin  n(j)'  J  rr'  dr'  d(f 


Llx  P,  J  (  2  ■) 

r^--n^ff^cos^'Vdr'^' 

oo  or 


(76) 


(77a) 


(77b) 


The  components,  Jr(r',tf>')  and  Jr(r',0')  are 
confined  in  the  horizontal  plane  to  a  small  area  bounded 
by  a  circle  with  radius,  r'max  =L«  l/kl.  In  this  case  we 
use  the  expansion 


/= CO  /=  CO 

J  =  'ZPM’)s\nl<jj  +  'ZQl(r')  coslfl 

1=0  1= 0 


(78a) 


/=00  .. - 

7/r\^')=  (/-')  sin/^'  +  ^(/-^cos/^'  (78b) 

/=0 


/= oo 

Si 

/=0 


/’VA,(rlr')./(r')dV'^Afaa(-A1)/M'oa(/»2)./(r')rfV' 

V  V 

(72) 


( h2)»J(r')  =  ^{ih2z2 ) 

nj  {Zr')  .  dJtl(Zr') 

- - sin  n(f)  ar - - - cos n(j)  a 

r'  r  dr' 


'VA+Jfy] 

(73) 


The  terms  denoted  by  J„(Zr')  again  refer  to  Bessel 
functions  of  integer  order  “«  ”,  and  are  not  to  be  confused 
with  the  components  of  current  J  r  and  J  , . 


The  radial  functions  of  Equation  (78)  are  analytic  with  no 
singularities.  Using  the  formulas 


2  n  In 

J  sin2  n(f>'d(p'=  J  cos2  rujZ  d(p'  =n 


Jsin  nt/>'cosm<ft'd(p]=0 
o 

we  get 

T %(Z)  =  -x]nJn(Zr')Pn(r')dr' 


(79a) 


(79b) 


(80a) 


Integrating  Equation  (73)  over  dV'  gives 


SM'enZ  (h2  W(?'  W  =  exp  (ih2z2  )rA 


r  =  r(1)  +  r(2) 

1  An  ~  1  An  ^  1  An 


rff  =-\nJ^Ar\m  n<f'J,.dV' 
v  r 


(74) 

(75a) 

(75b) 


t%W  =  -n)?±j^LSn{r'ydr' 

o  or 


rAla)=ca)+ca) 


(80b) 


(80c) 


Inserting  the  foregoing  results  into  Equation  (71a) 


gives 


Ea„ (?)  =  l^2,,  S<>)  cM e„,(-ht)exp(ih2z2)V AndZ  (81) 

0  /+7TAu1 


6 


2 


c 

h2 


/ij  +  h2 


(82) 


M  mA(-K)zx\>(ih2z2) 


exp  (;'/!,  z,  +ih2z2 ) 


nJJ.hr) 


sin  ntf>cir  + 


dr 


cos  nficip 


(83) 


EAn(r)  =  ~ 


cojuJ2-80)  _  ojfiJ2-SJ 


4/r 


-  K„a,  - 


4=1 

0 


^(/Zj  +  /z2 ) 

I  /\„  siting 


exp^'/ijZ,  +ih2z2 ) 


4n 
nJJAr) 


dA 


Lna t  (84) 


(85) 


Ln  =  J 


expp'/ijZ, 


/K/i,  +  4) 

F,„  cos  nii 


dr 


VIA 


(86) 


Following  the  procedure  leading  to  Equations  (84)  to 
(86),  we  readily  compute  the  remaining  terms  of  Equation 
(70)  and  thereby  compute  the  aboveground  field  to  be 
detected.  In  the  next  section  we  render  a  sample 
calculation. 


5.  COMPUTATIONAL  CONSIDERATIONS  AND 
SAMPLE  CALCULATION 


In  this  section  we  examine  the  challenges  in 
computing  the  integrals  in  Equations  (84)  to  (86),  and 
then  render  a  sample  calculation.  We  first  cast  Equations 
(84)  to  (86)  in  the  form 


EM(r)  =  - 


co/jJ2  -  S0)nsinn<fi  —  _ 


2  nr 


K.,a,.  - 


co/u0 (2  -S0) cos n<j>  —  _ 
2  n  "  * 


(87) 


field,  as  given  by  Equation  (70),  must  satisfy  the  radiation 
conditions  at  z{  — >  °°  and  z2  — »  oo  .  That  is,  the 
electromagnetic  field  must  go  to  zero  at  these  limits. 
Concern  about  this  arises  from  the  two  branch  points 

connected  with  h{  =  -J k*  -  J  and  h2  =  ^ k2  -  A 2  because 

of  the  ambiguity  of  taking  the  +  or  -  square  root.  These 
branch  points  each  then  lead  to  two  branch  cuts,  and  each 
branch  cut  leads  to  two  Riemann  surfaces.  Thus,  we  have 
a  total  of  four  Riemann  surfaces  to  consider.  At  every  step 
of  the  numerical  calculation  the  radiation  condition  must 
be  satisfied. 

We  note  in  passing  that  even  though  there  is  a  “  A~‘ 
factor  in  both  integrands  of  Equations  (85)  and  (86),  this 
does  not  lead  to  a  singularity  or  a  pole.  In  the  limit  of 
A—> Owe  have  Jn (Ar)  -» X  ,  so  that  the  only  term  that 

could  yield  a'  in  this  domain  is  J 0(Ar)  .  But  as  we  see 
from  Equation  (85),  the  factor  “n”  ensures  that  this  whole 
term  vanishes  for  n  =  0 .  We  also  have 


dJJAr)  xdJJg) 
dr  dg 


(90) 


where  g  =  Ar .  Since 


^ ^  is  always  finite,  the  A  term 


from  Equation  (90)  cancels  A  1  and  the  apparent  pole  or 
singularity  disappears. 


We  have  not  yet  attempted  a  numerical  evaluation  of 
Equations  (88)  and  (89).  Such  a  computation  will  be 
considered  in  the  event  that  a  near-term  analytical 
approach  based  on  the  “Method  of  Steepest  Descent” 
(MSD)  is  not  fruitful.  The  mathematical  structure  of  the 
integrands  of  Equations  (88)  and  (89)  is  consistent  with 
that  used  in  MSD  (Chew,  1990)  and  in  the  Ohio  State 
mathematics  lecture  series.  Let  S  represent  either  Kn  or 
Ln  since  they  are  equivalent  in  structure,  Q(  A)  represent 


either 


A(ht  +h2) 


rA,  dJJAr) 
A{hi  +  h2 )  dr 


and 


f(A)  =  ihlzl  +ih2z2. 


(91) 


K„=\ 


JJ^)ea„ 

A(ht  +h2 ) 


expO'/ijZj  +ih2z2 ) 


dA 


(88) 


4=1 

0 


dJJAr) 


A(hl  +hj  dr 


expO'/ijZ,  +ih2z2) 


dA 


(89) 


The  MSD  form  of  Equations  (88)  and  (89)  is 

oo 

S  =  fQ(A)expf(A)dA  (92) 

o 

Branch  points  from  h]  (A)  and  hJA)  require  integration 
in  the  complex  plane  using  A=x  +  iy .  Asymptotic 
approximations  to  Equation  (92)  require  a  number  of 
conditions.  First,  we  must  have  the  situation  where 
Q(A)  is  a  slowly  varying  function  of  A  in  the  region  where 


The  foregoing  equations  are  in  a  suitable  form  for 
numerical  integration  provided  we  are  aware  of  the 
guiding  principle  that  the  sum  of  all  constituents  of  the 
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exp  /(A)  is  the  major  contributor  to  the  integral.  Next,  we 
require  that  /(A)  =  u{x,  y)  +  iv(x,  y)  be  analytic  and  the 
functions  u(x,  y)  and  v(x,  y)  satisfy  the  Cauchy-Riemann 
equations.  Assume  that  a  saddle-point  extremum  of 
/(A)  exists  at  a  point  A0  determined  by  solving  the 
equivalent  equations 


far 

=  0  = 

k 

CONCLUSIONS 

Using  the  dyadic  Green’s  function  approach  this  study 
provides  a  rigorous  methodology  for  determining  a 
detectible  aboveground  electromagnetic  field  for  a  class 
of  buried  target  currents.  Numerical  computational 
techniques,  as  well  as  totally  analytical  techniques  based 
on  the  asymptotic  theory  of  integration,  are  examined. 
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Using  Equation  (93)  we  have  the  following  result  in 
the  neighborhood  around  A0 
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